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Abstract. We discuss the thermodynamic and dynamical properties of a variable dark energy 
model with density scaling as oc (1 + z)'", z being the redshift. These models lead to the 
creation/disruption of matter and radiation, which affect the cosmic evolution of both matter 
and radiation components in the Universe. In particular, we have studied the temperature- 
redshift relation of radiation, which has been constrained using a recent collection of cosmic 
(— I _ microwave background (CMB) temperature measurements up to 2 ~ 3. We find that, within 

^ ^ the uncertainties, the model is indistinguishable from a cosmological constant which does not 

I ■ exchange any particles with other components. Future observations, in particular measurements 

2 , of CMB temperature at large redshift, will allow to give firmer bounds on the effective equation 

■ of state parameter lUcff for such types of dark energy models. 



1. Introduction 



>: 

oo. 

Tjlj- . The current standard cosmology model relies on the existence of two unknown dark components, 

CN ■ the so called "dark matter" (DM) and "dark energy" (DE) , which amount to ~ 25% and ~ 70% 

of the total energy budget in the Universe, respectively. According to several observations, the 



I Universe is spatially flat and in an accelerated phase of its expansion [U [H [3l HI [5] . DE, described 



as a cosmological constant A in its simplest form, is modelled by a fluid with a negative pressure, 
which is a fundamental ingredient to explain the actual accelerated expansion of the Universe. 

Several models have been proposed to explain DE [H [3 El [9l [lOl [TTl [121 US] ■ An alternative 
consists to consider a phenomenological variable DE density with continuous creation/disruption 
/S '. of photons [a da [111 da HBl US] or matter [20l [2l]. The DE might decay /grow slowly 

■ in the course of the cosmic evolution and thus provide the source/sink term for matter and 

radiation. Different such models have been discussed and strong constraints come from very 
accurate measurements of the cosmic microwave background (CMB) radiation and other typical 
cosmological probe. 

CMB radiation is the best evidence for an expanding Universe starting from an initial high 
density state. Within the Friedmann-Robertson- Walker (FRW) models of the Universe the 
radiation, after decoupling, expands adiabatically and scales as (l + z), z being the redshift [22]. 
If we assume that each component is not conserved, contrarily to the standard scenario, then 
depending on the decay mechanism of the DE, the created photons could lead to distortions in 
the Planck spectrum of the CMB, and change the evolution of its temperature. The chance to 
appreciate the deviation from the standard temperature evolution is given by the increasingly 
number of recent works collecting observations of CMB temperature both at low [2'6\ [21] and 
higher redshifts 



Following the theoretical lines of |14 ^ 115 1 \17\. in Jetzer et al. [181 119j we have discussed a 
variable DE model A.{z) oc (1 + z)™' decaying into photons and DM particles. In particular we 
studied thermodynamical aspects in the case of a continuous photon creation, which implies a 
modified temperature redshift relation for the CMB. We have tested the predicted temperature 
evolution of radiation with some recent data on the CMB at higher redshift from both Sunyaev- 
Zel'dovich (SZ) effect and high-redshift QSO absorption lines as well as with an updated 
collection of data from different kind of observations, like distance moduli of Supernovae la, 
observations of the CMB anisotropy and the large-scale structure, together with observational 
Hubble parameter estimations [TO]. A similar approach, however without considering photon 
creation and thus a modification in the CMB temperature evolution, has been discussed by Ma 



2. Theoretical framev^ork 

We assume a cosmological framework based on the usual Robertson- Walker (RW) metric element 
|22j and that the Universe contains three different components: a) a matter (both baryons and 
DM) fluid, with equation of state pm = (since pm « Pm), b) a generalised fluid with pressure 
p-^/ = (7 — where 7 is a free parameter, which is set to 4/3 for a properly said radiation 

fluid, and c) a DE, x component, with pressure px and density p^- The equation of state for 
the X component could assume a very general expression, but we limit ourselves to consider the 
simple linear relation px = WxPx- We will set any 'bare' cosmological constant Aq equal to 
[18] . With these components we get for the Einstein field equations [22] 

k 

SirGptot = 3-^ + 3-^ , (1) 

S,Gp„ = -2| - I - ^ , (2) 

where ptot = P'y + Px and ptot = Pm + P7 + are the total pressure and density, R is the scale 
factor, fc = 0, ±1 is the curvature parameter and a dot means time derivative. Furthermore, we 
will assume that there is no curvature, thus A; = [3] . 

In the following we will adopt Wx = —1, however since we are assuming that the vacuum 
decays into radiation and massive particles, the effective equation of state Weff, which is the 
measured quantity, can differ from —1. 

Following [181 HS] (but see also [IH [16] ) we set here the energy conservation equation for the 
different fluids. The fluid as whole, verifles the Bianchi identity, which means that energy and 
momentum are locally conserved: V^T^'^ = 0, with T^'^ the stress-energy tensor 

T'^" = {ptot + Ptot^^W" - ptotg"" (3) 

and being the four velocity. After easy calculations, this conservation law reads 



ptot + ^{ptot+Ptot)H = , (4) 

where H = R/R is the Hubble parameter. In the standard approach, each component is 
conserved, thus Eq. ^ holds for all the fluid components, but here we will suppose that 
each component will exchange energy with each other. In particular, for matter, radiation and 
DE we impose the following relations (see also [26]) 

Pm + ^PmH = {I - e) Cx , (5) 



Py + 3'y pyH = € Cx , 



(6) 



Px + "iiPx + Px)ii = -Cx , 



(7) 



where we assume that both the matter and the radiation fluids exchange energy with the DE 
as parameterized by Cx and e. In particular, Cx depends on the DE and, indeed, acts as 
a source/sink term for the fluids energy. Evidently, if no interaction between the different 
components exists, then Cx is null and the standard picture is recovered. Moreover, if e = 
(e = 1), then the DE exchanges all the energy with matter (radiation). Cx can describe different 
physical situations such as, for instance, a thermogravitational quantum creation theory [TS] or 
a quintessence scalar field cosmology [7]. As we will discuss later on, e has to be very small 
(i.e., e <C 1), otherwise the radiation density would become much too big, contrary to present 
values. As an order of magnitude estimate we expect e ~ Ei+fi for which indeed e ^ 1, since 

+ P-y < Pm- 

Adopting as mentioned above the relation px = —px and defining px = A(t)/(87rG), from Eq. 
d?]) we obtain 

^- <«) 

We assume a power law model for the A{t) function, A(t) = B{R/Rq) or equivalently in 
terms of redshift, A(z) = B{\ + a;)™", where S is a constant. The value of this constant is 
B = ml{l-^mo), which can be found using Eq. ([1]) at the present epoch, assuming that today 
P7 ^ PmiPx- Thus, the density evolution for the x component is given by 

Px/Pcr^t = {l-^n.0){l+Z^, (9) 

where we have defined as pcrit = |^ the present critical density of the Universe. If m is 
positive, then the DE slowly decreases as a function of the cosmic time, whereas if m is negative 
the inverse process happens. 

From Eqs. ^ and ([6]) we derive the evolution laws for matter and radiation, 

Pm./ Pcrit = ^^mo(l + zf (10) 

m — 6 



P-t/ Pcrit = ^^70(1 + zf^ (11) 

_^m(l-a^ 3 ^ 

m — ^7 

where fimo and O^o are the matter and radiation energy densities at z = 0, respectively. 

Since we are interested in the evolution of the radiation temperature, it is useful to discuss 
the extreme case when only photons enter in the process (i.e. e = 1). Then for the matter 
density in Eq. (jlOp the usual evolution oc (1 + z)^ holds, while the radiation besides the usual 
term 0^o(l + z)'^^ has also a perturbative term depending on m. Since today El ^^70 ~ 5 x 10~^ 
it turns out that m has to be extremely small ^ 10~^. Therefore, unless m is extremely small 
or vanishing, DE has to decay mainly in matter with possibly some photons as well. Thus the 
condition e ^ 1 has to hold. 

^ The present radiation density is the only cosmological parameter accurately measured. The radiation density is 
dominated by the energy in the cosmic microwave background (CMB), and the COBE satellite FIRAS experiment 
determined its temperature to be T = 2.725 ± 0.001 K [27], corresponding to Q.^o ~ 5 x 10~^. 



2.1. Hubble and deceleration parameter 

Due to the very small value of Jl-yo it follows that the evolution of the Universe is essentially 
driven by the DE and DM, therefore, from Eq. ([1]) the following law for the Hubble parameter 
holds 



H[z) ~ -^{Pm + Px) 



3 
Ho 



3(1 - n^o) ^ ^ (30mO - m) ^ ^^3' 



m 



m 



1/2 



which is obviously the same expression found in Ma |21j . 
Recasting Eq. ([7]), it is possible to write 

(J 

p^ + 3H{p^ + Px + ^) = , 
which shows that the term Cx contributes to an effective pressure 

PcS 



Cx Cx 



Therefore, we get an equivalent effective DE equation of state Wes [21 

PcS m 



■WeS 



1. 



Px 



(12) 



(13) 



(14) 



(15) 



If m > then we have Wcs > — 1, i-e. our model is quintessence- like [TJ [T2], while we have a 
phantom-like model when m is negative and < — 1. Another interesting quantity is the 
deceleration parameter, which can be written as 



(1 + zf^m - Snmo) + 3(m - 2)(1 + z)'"(l)^o - 1) 
2(1 + z)3(m - Sn^o) + 6(1 + z)"'{n^o - 1) 



(16) 



Imposing that q{z) = 0, we can determine the transition redshift, i.e. the redshift at which the 
Universe changed from a deceleration to an acceleration phase, which is given by 



ZT 



3(2-m)(l-l)^o) 

30jnO ~ ITT- 



1 

3 — m 



1 . 



(17) 



Prom this result we see that the larger m is, the earlier the Universe changes from deceleration 
to acceleration. 



2.2. Thermodynamical aspects and CMB temperature evolution 

We follow here the approach outlined in Lima [H], where he defines a current as A'^" = nn" 
with n being the particle number density of the photons or of the DM particles. Indeed, there 
is a current for each of these components. Due to the decaying vacuum the current satisfies the 
following balance equation (one for each component) 



hi + ZriiH = ipi , 



(18) 



with i = 7 or DAI (7 for the photons and DM for the dark matter) and "01 is the corresponding 
particle source. For decaying vacuum models "07 + V'dm is positive and related to the rate of 
change of px- We can also define an entropy current of the form 



E 



(19) 



where C7j is the specific entropy per particle (photons, DM and in principle also DE). If the DE 
Px is constant the above entropy current is conserved. The existence of a non equilibrium decay 
process of the vacuum implies S^^ > 0, thus an increase of the entropy as a consequence of the 
second law of thermodynamics. In principle the second law should be applied to the system 
as a whole, thus including the vacuum component [26]. Assuming that the vacuum is like a 
condensate with zero chemical potential p^ac it follows from Euler's relation 



Pv 



Px +Px 



n 



(20) 



provided Wx = —1, that a^ac = and thus its contribution to the entropy current vanishes. 
Given our assumptions the vacuum plays the role of a condensate carrying no entropy. 

For instance, for quintessence models in the limit where the scalar field does not depend on 
time, and thus its time derivative vanishes, one gets Wx = In the later stages of the Universe 
the time dependence is possibly very weak so that Wx = holds up to small corrections. 

The equation for the particle number density of radiation component is given by Eq. (jlSp 
with i = J. Using Gibbs law and well-known thermodynamic identities, following the derivation 
given in the paper by Lima et al. [T6j, one gets (see also [HI [18]) 



Pj + Py 



^7 



(21) 



To get a black-body spectrum the second term in brackets in Eq. ()2ip has to vanish, thus 



Thus, Eq. ([21]) becomes 



eCx = — [p-y + P7] . 



T _ fdp. 



With ^1^^ = (7 ~ 1) one obtains 



T 

T 



(7-1 



(22) 



(23) 



(24) 



Using the equation for the particle number conservation Eq. ([T8]) into Eq. ([2^ leads to 



T 

^ = (7-1) 



■07 



(25) 



With Eqs. ([22D and ([25D we get 



^ = (7 



eA 



8TrG{p^ + p^) 



3H 



(26) 



Now, following the previous discussion on e and aiming to be very general, we set e = Billii^ 
where e is a new parameter and insert it into Eq. ()26p . Taking the sum of Eqs. ([T]) and ([2]) we 
find 

87rG(ptot + ptot) ^ STTGpm = 2j^-2- = -2H . (27) 



R 



Finally, we obtain the expression 



T 
T 



(7-1) 



which we can integrate 



dt=ij- 1) 



Ae 








-3H 


1 


2H 








" Ae 






?,H 


hi 


2H 



(28) 



dt 



(29) 



where to denotes the present time and ti some far instant in the past. Indeed, if A vanishes and 
7 = 4/3 one gets the usual dependence T{t) = ^^^ji^j^^ ^ for a radiation fluid. To carry out the 
integration of the first term on the right hand side it is useful to perform a change of variable 
from t to z and accordingly ^ = jf^Y+I)' This way we get (with zi corresponding to the time 
ti and zq = corresponding to to present time) 



In 



T{z = 0) 

nzi) 



+ 3(7- l)/n 



(7-1) 



R{z = 0) 

R{zi) 
A'e 



(30) 



H'H{1 + z] 



dz 



where ' denotes derivative with respect to z. 

As next we insert H{z) and its derivative as taken from Eq. (jl2p into Eq. ()30p and integrate 
it, to get (setting zi = z) 



T{z) 



Rn 



R{z] 



3(7-1) 



exp 



B{1 - lYe 



-A 



mO 



(31) 



where 



A = ln{{m - Sn^o) + "i(l + z)"'-'\n^o - 1)) 

-ln{{m - 3)Qrno)] ■ 



(32) 



We can also write Eq. ([31]) as 

r(z) = ro(i + z)3(T-i) 

(m - 3n,no) + m{l + z)"^-3(0 



mO 



(m - 3)0 



mO 



?(7-l) 



(33) 



We inserted in the exponent of Eq. ()3ip the explicit form of B, thus getting as exponent in 
the above Eq. e(7 — 1). Hereafter, we will set e = 1. Clearly e and m are not independent, 
we checked using the temperature redshift data that if e is bigger (~ 10 or more), then m has 
to be extremely small consistently with what mentioned above (as it would lead to a too high 
production of photons in the DE decay). On the other hand, if e gets smaller (e.g., ~ 0.1) m 
gets bigger (~ 0.2) and accordingly w^fj, moreover m would be poorly constrained, since the 



uncertainties would then be very high. But from the other data (without the the temperature 
ones) there are aheady stringent hmits on m and thus this way one could get lower limits on e, 
under the assumption that DE decays also in photons. 

Notice that for z = we have T(0) = Tq, whereas for m = the expression in the parenthesis 
is equal to 1 and thus T{z) = Tq{1 + z)^^'^~^\ which for the canonical value of 7 = 4/3 reduces 
to the standard expression. 

3. Results 

To test the temperature evolution for the radiation component, we rely on the CMB temperatures 
derived from the absorption lines of high redshift systems and the ones from SZ effect in 
clusters of galaxies (we will collectively quote as Tcmb, hereafter). At high redshift the CMB 
temperature is recovered from the excitation of interstellar atomic or molecular species that 
have transition energies in the sub-millimetre range and can be excited by CMB photons. When 
the relative population of the different energy levels are in radiative equilibrium with the CMB 
radiation, the excitation temperature of the species equals that of the black-body radiation at 
that redshift, providing one of the best tools for determining the black-body temperature of 
the CMB in the distant Universe [28l[29l[30l[3ll[32l[33l[Ml[35]. We used a sample of 9 QSO 
absorption measurements |19^ [25] . In summary we have 4 data points from the analysis of the 
fine structure of atomic carbon (AC) and 5 measurements based on the rotational excitation of 
CO molecules (CO) [25]. 

At lower redshift we use the measurements from the SZ effect. Thus, spectral measurements 
of galaxy clusters at different frequency bands yield independent intensity ratios for each cluster. 
The combinations of these measured ratios permit to extract the cosmic microwave background 
radiation (see Fabbri et al. [36] )• We relied on the data compilation in Luzzi et al. [23], which 
have analyzed the results of multifrequency SZ measurements toward several clusters from 5 
telescopes (BIMA, OVRO, SUZI II, SCUBA and MITO). 

We derive the observed Tqmb from the theoretical expression Tffi, which we have derived in 
Eq. ([55]) . We set Tq = 2.725 K, which is quite well determined in the literature [33, and the 
matter density r^mo = 0.273 to the value inferred in Komatsu et al. [SS]- If we take 7 = 4/3, 
then we find m = 0.03t[!;J]^ |19j . which is lower than the estimated value of m = 0.09 it 0.10 in 
|18j . but fully consistent within uncertainties, and also pretty consistent with m = 0. In Fig. [1] 
the temperature measurements (together with some upper limits) are shown, and our best fitted 
result is plotted and compared with the m = —0.09 result in Ma |21J. The value we have found 
corresponds to an effective equation of state w^s = —0.99 ±0.03, consistent with Wcs = — 1> and 
the transition redshift is zt = 0.78 it 0.08. In order to check the impact of redshift distribution 
we separate the data in two redshift bins with z < 0.6 (SZ data only) and z > 0.6 (QSO 
absorption lines only), finding the best fitted values m = 0.12~^q'\^ and -Omtoll, respectively. 
Although the uncertainties are very high, it seems to emerge a mild trend with lower redshift 
data preferring a DE decaying into matter and radiation, while data at z > 0.6 point to an 
opposite behavior. These results could be interpreted in a different way, in fact, the differences 
found could be due not to the different redshift coverage, but to some particular biases in the 
two kind of observations, SZ vs QSO absorption lines. An indication of this suggestion comes if 
we divide the sample in three subsamples: 1) SZ data, 2) the data from the analysis of the fine 
structure of atomic carbon (AC), and 3) the measurements based on the rotational excitation of 
CO molecules in |25j . If we fit the model to the combined SZ+AC and SZ+CO samples we find 
m = 0.11~^q'W and m = 0.03l^g'j'g, respectively. Because of the larger measurement errors, the 
AC data affect very little the results when only SZ measurements are adopted, and the result 
for the second sample shows that SZ and CO data mainly constrain m, and AC simply gives a 
tiny reduction on the errors. 

Adopting a constant value for the ratio ip^/Sn^H = (3 Lima et al. |16j have found the simple 
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Figure 1. Cosmic microwave background temperature as a function of the redshift. The black 
points with bars are the full collection of measurements from Luzzi et al. [23] and Noterdaeme et 
al. |25] . The gray arrows represent the upper limits derived from the analysis of atomic carbon 
(see [25] for details). The black line is the best fit result (m = 0.03), while the gray region is the 
la uncertainty. The red line is the best fit recovered from Ma [21] (see the legend). The inset 
panel show a magnified vision of the higher redshift region of the plot. 

relation, 

T{z)=Toil + z)'-^ . (34) 

If we fit all the sample we find /? = —0.002 it 0.03, while for the low and high redshift subsamples 
we have /3 = 0.06 it 0.08 and —0.01 ±0.03, respectively. These results are qualitatively consistent 
with what found in Noterdaeme et al. [25], and points to a similar trend as the one discussed 
above. 

If 7 is left free to change, for the best fitted value it turns out that 7 > 4/3 and m is 

_i_n r\'i _i_n 

systematically more positive. We find 7 = 1.35_q q3 and m = 0.25_q'^7, which corresponds to 
an effective equation of state = —0.92 it 0.07 and the transition redshift is zt = 1.1 ± 0.6. 
When the two subsamples are adopted, wide confidence contours are found, particularly for the 
z > 0.6 sample, for which the contours at very low m are not closed. We obtain 7 = 1.3]!^q^ and 
1.261q'q}, while m = O.S^Qg and 0.6^5 qj respectively for the low and high-z samples. 

We notice that if DE does not decay into radiation (corresponding to e = e = and thus m 
is no longer constrained) then the CMB temperature will scale in the standard way. Clearly, 
this would imply that if DE decays, this has to be into DM only. On the other hand a deviation 
of the CMB temperature from the standard scaling could be interpreted as DE decaying also 
into radiation. In which case with Eq. ()33p one can determine either m or e and thus get some 
insights on the decay mode of DE into radiation. Future data on the CMB temperature will 
allow to shed light on this important issue. 

In |19j we have combined the CMB temperature data with a set of other independent 
measurements such as Supernovae la distance moduli, CMB anysotropy, baryon acoustic 
oscillation and observational data for the Hubble parameter. When combining all the data 
sets we found values for m, and thus w^s consistent with the ones obtained when using only 
CMB temperature data. 

C. Tortora was supported by the Swiss National Science Foundation. 
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